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Abstract 

Dirac equation written on the boundary of the Nutku helicoid space con- 
sists of a system of ordinary differential equations. We tried to analyze 
this system and we found that it has a higher singularity than those of 
the Heun's equations which give the solutions of the Dirac equation in the 
bulk. We also lose an independent integral of motion on the boundary. 
This facts explain why we could not find the solution of the system on the 
boundary in terms of known functions. We make the stability analysis of 
the helicoid and catenoid cases and end up with an appendix which gives 
a new example where one encounters a form of the Heun equation. 
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1 Introduction 



Although one usually needs only different forms of the hypergeometric equa- 
tion or its confluent forms to describe many different phenomena in theoretical 
physics, functions with higher singularity structure are seen more and more in 
the literature [1-12]. A common form is the Heun function [T3], which is studied 
extensively in the books by Ronveaux and Slavyanov et al [T3] [T3] , the seminal 
book by Ince Hi as well as in several articles [17-20]. Although for linear equa- 
tions both Ince and Slavyanov et al end their singularity analysis with Heun 
type functions, sometimes equations with even more singularities are needed for 
relatively simple situations, which lack some symmetries. 

As an example of such a case, here we study the singularity structure of the 
Dirac equations, written in the background of the Nutku helicoid solutionpTj. 
restricted to the boundary of the helicoid. This metric was formerly studied 
by Lorenz-Petzold One can study the scalar field in this background and 
obtain the propagator in a closed form [53], thanks to four integrals of mo- 
tion allowed by the metric. One can also write the Dirac equation and obtain 
the solutions in terms of Mathieu functions [24-26]. One needs to study the 
eigenvalue problem on the boundary to impose the boundary conditions in this 
problem. The similar problem in the bulk, using partial differential equations, 
can be solved in terms of known functions. On the boundary, we get a system 
of ordinary differential equations. At first glance, this system seems to be easier 
to analyze. When we investigate the system further, we find that this system 
has higher form of singularities and fewer integrals of motion. These turn out 
to be the reasons why we cannot express the solution in terms of the functions 
cited in luce's book[16j. 

In this note, we comment on the symmetries of the new problem and study the 
singularity structure of the new system. We show that one gets an equation of 
higher form of singularity. Since we do not have closed solutions, we then use 
stability analysis to see if this system describes a stable system. To our surprise 
we find that although the answer is not affirmative for the helicoid case, we get 
a limit cycle for the related catenoid solution. 

Below we summarize our results. In an appendix, we show a new example where 
one encounters a form of the Heun equation. This is the solution of the laplacian 
in the background of the Eguchi- Hanson solution , trivially extended to five 
dimensions. 

2 Our Analysis 

The Nutku helicoid metric is given as 
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a2 



ds^ = [dr' + jr'' + a')de^ +{ I + — sin' ] dy^ 

2 / 2 \ 

-^sm29dydz+( l + ^cos'^ 9 j dz'^]. (1) 

where 0<r<oo,0<6'<27r,y and z are along the KiUing directions and will 
be taken to be periodic coordinates on a 2-torus [23J . This is an example of a 
multi-center metric. This metric reduces to the flat metric if we take a — 0. 

ds^ = dr^ + r'd0^ + dy^ + dz^ . (2) 

If we make the following transformation 

r = asinhx, (3) 

the metric is written as 

a? 

ds' = — sinh 2x{dx' + dO'^) 

- sin29dydz + (sinh^ x + cos^ 0)dz'^]. 

The solutions of the Dirac equation, written in the background of the Nutku 
helicoid metric, can be expressed as a special form of Heun functions [25]|26j. 
We could reduce the double confluent Heun function obtained for the radial 
equation to the Mathieu function with coordinate transformations. Mathieu 
function is a related but much more studied function with similar singularity 
structure. In the work cited above, [25] [26], we tried to get the solution of the 
little Dirac equation, the name used for the equation restricted to a boundary 
of the helicoid. We needed this solution also to be able to calculate the index of 
the differential operator. For a fixed value of the radial coordinate we had only 
a coupled system of ordinary differential equations, which, in general, should be 
much simpler to solve than the coupled system of partial differential equations 
obtained for the full Dirac equation, written for the bulk. We were successful to 
obtain the solution in this latter case in terms of Mathieu functions. We were 
not able to identify the solutions for the little Dirac equation , though. 
The first thing we check is whether we lose any of the three Killing vectors and 
one Killing tensor. For the metric in question, one has three integrals of motion, 
|23| namely Py,Pz and g^^p^Pu — /i^, namely 

^("f^) "^"^ (pI+pI) sinh^a;-^^ sinh2a;+a2 (cos6'py + sine'p^)^ = 0, 

(5) 
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and an extra integral of motion, the Killing tensor [23] , 



K ~ —pg — a {cosOpy + sinOpz) (6) 

which gives us the angular equation for a fixed value of the constant A . See 
eq.s (46) and (47) of [23 . When one restricts the solution to a fixed value of 
the radial coordinate, the value of g^'^PfiPi, is not an independent constant of 
motion from the other two Killing vectors and the Killing tensor. 
Trying to see the effect of this result on our problem, we investigate the singu- 
larity structure of the equation we get for the little Dirac operator. Here we 
study the simplest case, where the eigenvalue A is equal to zero, since the an- 
swer to our problem is already apparent here. For this case, instead of getting 
a system of four coupled equations we get coupling only between two of them 
at a time. When we analyze the system by reducing them to a second order 
equation for a single dependent variable, we find an operator with two irregular 
and one regular singularities, which is one more than allowed for the equations 
considered among the Heun functions. The double confluent Hcun function, the 
solution one obtains for the full Dirac equation, has two irregular singularities, 
missing the extra regular singularity of the case studied here. 
To make our discussion concrete we explicitly perform the calculation in the 
next section. 

2.1 Singularities 

The Dirac equation written in the background of the Nutku helicoids metric is 
written as 

{dx+ide)'^^ +iak[cos{e-(l) + ix)]'S 4 = 0, (7) 

{d^-ide)-^4 -mfc[cos(6l-0-ia;)]*3 = 0, (8) 

i-d,+id8)h +iak[cos{e-cj) + ix)]f2=0, (9) 

{~dx-ide)f2 -iak[cos{6~-(f)~ix)]fi^Q. (10) 

These equations have simple solutions (23] which can also be expanded in terms 
of products of radial and angular Mathieu functions [2l][25]. Problem arises 
when these solutions are restricted to boundary [26]. 

To impose these boundary conditions we need to write the little Dirac equation, 
the Dirac equation restricted to the boundary, where the variable x takes a fixed 
value xi). We choose to write the equations in the form, 

^{j^*3 +ikacos{6~(j) + ixQ)^4} = \h, (11) 



a 



^{-*^*4 - iakcos{6 ixo)^^} - A/2, (12) 



a 



— -^akcos{e~cj, + lXo)f2} = X^^, (13) 

a da 
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— {*:^/2 +iakcos{0-^-ixo)fi} = X^4- (14) 
a au 

Here A is the eigenvalue of the httle Dirac equation. We take A = as the 
simplest case. The transformation 

e = e-(j)~ixQ (15) 
can be used. Then we solve /i in the latter two equations in terms of f2- 

{ak) 

(16) 



-j^/2-tane — /2+^^[cos(2e)cosh(2a;o)-isin(2e)sinh(2a:o)+cosh(2xo)]/2 = 



de^' de" 2 

When we make the transformation 

u-e^'®, (17) 

the equation reads, 

{4(u+l)M[u^ + ^]-2iM(?/-l)^+.^(ii+l)[ue-2-o + ie2-o+cosh(2xo)]}/2 = 0. 
du^ du du 2 u 

(18) 

This equation has irregular singularities at m = and oo and a regular singularity 

oo 

at li = — 1. If we try a solution in the form ^ a„M" around the irregular 

n— — oo 

singularity m = we end up with a four-term recursion relation as 

a„_i[4(n2 - 2n + 1) - 2^(n - 1) + M!(^e-2-o + Ig^.o)] 
+a„[4n2 + 2^n + ^(^e"^-" + \e^-^')] + (19) 

As it is known, in the Heun equation case, this kind of series solution gives a 
three-term relation [29] , 

If we search for a solution of the Thome type we may try a solution of the form 

A 

f2 = e^g{u). This form docs not allow us to get a Taylor series expansion 
around the irregular point u = [5] [30] . 

If we try a series solution around the regular singularity at u = — 1 as 

C30 

^ a„(u -I- 1)"+" we find a relation between five consecutive coefficients for the 

solution. Therefore, we may conclude that the solution of this equation cannot 
be written in terms of Heun functions or simplier special functions. 
To check this further, we first set the coefficient of ^ term in equation [T51 equal 
to zero to change our irregular singularity at zero to a regular one. Then we 
keep this term and discard the ue~'^^° term to reduce the singularity structure of 
infinity. In both cases one can check that the solution can be expressed in terms 
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of confluent Heun functions. This shows that reducing one of the singularities 

yields a Heun function. Thus, we conclude that the full equation [18] is not one of 

the better known equations in the literature, which are included in the computer 

packages like Maple, cited in the seminal book by Ince [TO] , 

To investigate the type of our equation we try to get a confluent form of a new 

equation, 

V {z) + ( + + )v (z) + Y( ^ y^^i = (20) 

z z + L z — a z'^(z — a} 

with regular singularities at 0, —1 and a and an irregular singularity at infinity. 
This equation differs from the generalized Heun equation [7]|17|: 

z z + 1 z — a z(z + l)[z — a) 

which also has regular singularities at 0, —1 and a and an irregular singularity 
at infinity. These two equations both have four-term recursion relations. They, 
however, have different singularity ranks according to the classification given in 
[31] . When we put a = in equation [20] we get 

/// X , /2 - Mo - M2 1 - Mu // X , (3q + I3iz + I32z'^ , . „ .„„x 

V \z) + ( + , -I )y (z) + 5 v{z) = 0. (22) 

z z+1 z-^ 

We get a singularity structure as a regular singularity at —1 and two irregular 
singularities at zero and infinity like the equation 1181 

Both equations [18] and [22] have four-term recursion relations when a Laurent 
power series solution is attempted. We may name this equation as the confluent 
form of the equation l20l It is in the same form as our original equation rewritten 
as 

+ n) + ^l/ + + + = 0. ,23) 

du^ u 2 u + L du 2 u 

Both of these equations have s-rank multisymbols {1, |, |} referring to the sin- 
gularities at {— 1,0,00} [31]. 

We could not obtain a confluent equation similar to the equation [18] from the 
generalized Heun equation l21l If we simply put a = in this equation we get the 
confluent Heun solution. We can obtain an equation with the same singularity 
structure as our equation only if we write the equation, 

"( A , ^1 ^ , 1 ^ ^1 , 1 " M2x u . . ^ + + I3iz + I32Z'^ 

y r+ )y{z)+— — r — v{z) = Q, (24) 

z z + 1 z~a zyz + l)[z — a) 

and make a approach zero. Then we end up with an equation having the same 
singularity structure as in equation 1181 

If we want to compare equation [TS] with equation]^ we have to form a confluent 
form of the latter equation. To coalesce the singularities at zero, we make a 
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detour and then use standart techiques [3^. We first translate the singularity 
at zero to a singularity at b, 



1 - Mo ^ 1 - A*i ^ 1-^2 
z ~ b z + I z ^ a 



-a)y'{z)- 



(3o + Piz + (32Z^ 
[z - b){z + l)(z - a) 



y(z) - (25) 



This equation has regular singularities at z — ~l,a,b and an irregular singu- 
larity at infinity. We make the transformation z = ^. Then, the equation [55 
becomes, 



y"{v) + (— 

1 - Mo 1 



(1 - ^o) - (1 - Ml) - (1 - M2) 



Ml 



V 

M2 



00 



+ 



)y'{v) 



{I - bv){l + v){l - av) 



y{v) = 



We set /io — ~Mi = l/b = 2/a — e and take the limit e 
back to the original variables using v — - to obtain 



(26) 



00. Then we transform 



y"{z) + ( 



3 - Ml 



Ml 



z{z + 1) 2' 



r)y'{z) 



z2(z+l) 



y(2) = (27) 



This equation is a confluent form of the equation [551 It has a regular singularity 
at —1 and two irregular singularities at zero and infinity and a four- way recursion 
relation when expanded around zero using a Laurent expansion. This equation 
has s-rank multisymbols {1,2,2} referring to singularities at {—1,0, 00} [31] . 
Here we get the rank-2 irregular singularities at zero and infinity only from the 
coefficient of the first derivative whereas in equation [TS| and in equation [2H the 
coefficient of the term without derivatives gives us these singularities. They 
also have different ranks. Even if we set a = in equation [57] then we get 
rank={l, 2, |} which is different from the rank of our original equation. Hence, 
our equation 1 181 may be a confluent form only of a variation of the equation [211 
like equation [24l 



2.2 Stability analysis for the little Dirac Equation 

The little Dirac equation is a system of linear differential equations with periodic 
coefficients. Then we can write the system as [33] [ 34 j : 

deij = Pi0)^P (28) 

Here P{9 + t) = P{9) and r is the period of the coefficients, (r 7^ 0). According 
to Bellman, if "0(0) = / , we can write 

m = Q(0)e«^- (29) 
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The matrix Q{9) is also periodic with the period r. Then we have 



Be„BT 



^{e + T) = Q{6 + T)e'"'e 



(30) 



■ip{0)e 



Bt 



We can define e = C and use 



C = VW"V(e + r) (31) 
to obtain C. The Jordan normal form of C, T being the transformation matrix, 



C = T 



\ 



Lr) 



gives us the B matrix: 



B = -\nL 

T 



(32) 



(33) 



The eigenvalues of B give us the characteristic roots. We will use these charac- 
teristic roots with different parameters in our stability analysis. The character- 
istic roots are given by aj, (i = 1..4). 

Table 1. The change in the characteristic roots with respect to parameters 

(helicoid case) 



a 


k 




A 


real parts of the 
characteristic roots ( x 27r) 


signs of the 
characteristic roots 


1 


1 


1 


1 


Q!i,2,3,4 = 9.25029 


+ H 


0.5 






11 


ai,2,3,4 = 4.03926 




0.8 


31 


11 


11 


ai,2,3,4 = 7.2163 




1.1 




11 


11 


Q!i,2,3,4 = 10.2542 




1.5 








"1,2,3,4 = 14.2215 




1 


0.5 


1 


1 


ai,2,3,4 = 5.8327 




J) 


0.8 






ai,2,3,4 = 7.75421 




» 


1.1 


11 


11 


"1,2,3,4 = 10.0319 






1.5 






0:1.2.3.4 = 13.2745 




1 


1 


0.5 


1 


ni.2.:!,l = (). 50200 








0.8 




"1,2,3,4 = 7.91586 






11 


1.1 


11 


"1,2,3,4 = 10.0632 






11 


1.5 


11 


"1,2,3,4 = 14.4648 




1 


1 


1 


0.5 


"1,2,3,4 = 8.30164 




5! 


11 


11 


0.8 


"1,2,3,4 = 8.81342 








11 


1.1 


"1,2,3,4 = 9.49249 






•>•) 


11 


1.5 


"1,2,3,4 = 10.5888 
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Our calculations indicate that /i and /2 solutions arc not stable (positive char- 
acteristic root) while, "^3 and "^4 solutions are stable (negative characteristic 
root). As it is seen in Table 1, when we keep all the other parameters constant 
and vary only a, the value of the roots arc influenced most, whereas the effect 
of the variation in the value of A changes the value of the roots least. We also 
find that when these parameters exceed unity in absolute value, we encounter 
inconsistencies in the numerical values. The separation between conseciitivc 
roots increase and some negative roots go to positive values for large values of 
the parameters. If we keep the values of the parameters in the range [—1, 1], we 
seem to have no such problems. 

The periodicity of the defined Q can be checked using numerical means. 
We use 

^{e + t) = Q{9 + T)e''^C (34) 
and for 6» = 0, ip{e) = Q((9)e^^ to give, 

i^{T)C-^ = Q{t) = Q(0) = V(0) = I. (35) 

We check numerically that this equation is satisfied; hence, Q is periodic. 

For the catenoid case, one replaces a with ia in the metric. The same stability 
procedure is performed for this case and we find the characteristic roots given 
in the Table 2. 

Table 2. The change in the characteristic roots with respect to parameters 

(catenoid case) 
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a (x i) 


k 


Xq 


A 


characteristic roots ( x 27r) 


signs of the 
characteristic roots 


1 


1 


1 


1 


ai 2 = -5.84543 x 10"** + 2.08426« 
0-3^4 = -6.40967 X 10-** + 0.72266i 


H h- 


0.5 




55 


55 


ai,2 = -2.40553 x lO"'^ - 1.19724i 
a3,4 = -2.49225 x lO"'^ - 2.3721i 


- H h 


0.8 


75 


55 


55 


ai.2 = -5.63046 x 10"^ - 0.233835i 
a3^4 = ~5. 78178 x 10"^ - 2.42696i 


- H h 


1.1 


15 


55 


55 


ai^2 = -7.82856 x 10" ' + 2.9982i 
a3^4 = -7.88358 x 10"^ - 0.186969i 


H h- 


1.5 


15 


55 


55 


ai,2 = -7.53539 x 10" ' - 0.526163i 
asA = -7.7255 x lO"'^ - 2.35978i 


- H h 


1 


0.5 


1 


1 


ai,2 = 3.49081 x 10"» + 0.661356i 
= 2.02794 x 10"** - 1.31726i 


H h 




0.8 


55 


55 


ai^2 = -4.04957 x 10"^ + 0.864554i 
asA = -4.22823 x 10""^ + 2.3638i 


H h- 




1.1 


55 


55 


ai,2 = -2.2577 x 10"^ + 2.7328i 
asA = -2.30929 x 10"'^ + 0.120058i 


+ - +- 




1.5 


55 


55 


ai.2 = -1.33916 X 10"" + 0.715427i 
a3,4 = -1.38904 x 10"^ - 2.75007i 


H h 


1 


1 


0.5 


1 


ai 2 = -4.68625 x 10"** - 0.232525i 
a3^4 = -5.34771 x 10"** - 2.02431i 


- H h 




55 


0.8 


55 


ai 2 = -4.69901 x lO" "" - 1.06977i 
^3^4 = -4.72986 X 10"^ - 2.11709i 


- H h 


55 


55 


1.1 


55 


ai.2 = -6.06425 x 10"'^ - 2.74153i 
a3,4 = -6.14188 X 10"^ - 0.143986i 


- H h 




55 


1.5 


55 


ai,2 = -4.26375 x 10"^ + 0.299098i 
a3^i = -4.4214 x lO""^ + 1.55822i 


H h- 


1 


1 


1 


0.5 


ai.2 = -4.67773 x 10"^ + 0.761149i 
a3^4 = -4.69627 x 10"'^ + 0.741915i 


H h- 


!! 


55 


55 


0.8 


ai,2 = -2.82256 x 10"^ - 1.5208i 
a3'4 = -2.86571 x 10"^ - 0.79892i 


- H h 


!! 


55 


55 


1.1 


ai,2 = -4.14479 x lO"'^ + 2.38047« 
a-i^^ = -4.1985 X 10"^ + 0.650337i 


H h- 


!5 




55 


1.5 


ai.2 = -3.33914 x lO""^ - 2.62503i 
a3^4 = -3.42067 x 10"'^ - 0.145628i 


- H h 



We see that the real parts of these roots are compatible with assigning to value 
zero within numerical errors. This corresponds to a limit cycle [35] [36]. 

3 Conclusion 

Here we performed a systematic analysis of the Dirac equation restricted to the 
boundary when it is written in the background of the Nutku helicoid solution 
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P5] We find that the resulting system of ordinary differential equations has a 
singularity which is higher than those of the Heun functions which are solutions 
for the bulk. We also lose an independent integral of motion. This fact explains 
why we could not obtain the solution of the system on the boundary in terms 
of well known functions. 

The stability analysis we performed shows that although this system is not 
stable, a related system, the catenoid solution is. We can, thus, give a meaning 
to its solutions, although we can not get explicit solutions for the little Dirac 
equation obtained from it too. 



4 Appendix: Scalar field in the background of 
the extended Eguchi-Hanson solution 

To go to five dimensions, we can add a time component to the Eguchi-Hanson 
metric [571 so that we have 

ds' = -de + -^dr^ + r^al + a^) +r'{l- ^)al (36) 

where 

ax = cosfrf^ — sin 9 sin ^dcj)) (37) 
ay = -{sinS^dO — sin 6 cos dip) (38) 
a^ ^^{-d£_- cos ed(l)). (39) 

This is a vacuum solution. 
If we take 

wc find the scalar equation as 



Hip{r, 9) = ( :^drr + ^ 9r + k^r^ + T + 

r r'' a — r 

, 8mncos0 - 4(771^ + n^) , , ^, ,^ , 

4906 + 4 cot 6*96 + —ri (41) 

sm~ 

If wc take (p{r,6) = f{r)g{9), the solution of the radial part is expressed in 
terms of confluent Heun (He) functions. 



/ (r) = {-a^ + r^) ^ " He [o, m, m, \ k^a' , i - i A - 1 k'a^ 
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+ {a^ + r^)-^ " (r^ - a^) ^ ffc (o, -m, m, \ k'a^ i - i A - J fc^a^ ~ 
The angular solution is in terms of hypergeometric solutions. 

g (0) = ^{^2-2 003^ (i cos (0) - ' " cos ^ - 1) 

[(2 cos (6*) + 2)^ 



2a2 
(42) 



12 2 " 2 ' 



X 2F1 {[-n + ^ VA + l + ^. - ^ VA + 1 + i], [1 - n - m], i cos {9) + i) 

(43) 

+ (2 cos (0) +2)^+5"+^™ 

X2fj([TO+ivA+l + i,m-^VA+l + i], [l + n + m],i cos(0) + i)]} 

If the variable transformation r = a\/ cosh a; is made, the solution can be ex- 
pressed as 



/ (x) = ^.^"^^^ {(sinh (a;))™"*"^ He ^0, m, m, ^ fc^a^, ~ ^ h^o?., ^ cosh (x) 

+ (2 cosh (x) + 2)^5 "+ 5 (2 cosh (x) - 2)^ ™+5 (44) 

2. 2 

We tried to express the equation for the radial part in terms of u = to 
see the singularity structure more clearly. Then the radial differential operator 

reads 

(f . . f I . 1\ d . ,2 2 / 1 . 1\ . 



dv? \u—l'^u)du'^^^ \u— 1'^ uj ^ M^(l — uY ^ 

This operator has two regular singularities at zero and one, and an irregular 
singularity at infinity, the singularity structure of the confluent Heun equation. 
This is different from the hypergeometric equation, which has regular singular- 
ities at zero, one and infinity. 
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